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Abstract 

The decay rate for a particle in a metastable cubic potential is investigated in the quantum 
regime by the Euclidean path integral method in semiclassical approximation. The imaginary 
time formalism allows one to monitor the system as a function of temperature. The family of 



O 

ri , classical paths, saddle points for the action, is derived in terms of Jacobian elliptic functions whose 

_. periodicity sets the energy-temperature correspondence. The period of the classical oscillations 

^ • varies monotonically with the energy up to the sphaleron, pointing to a smooth crossover from 



the quantum to the activated regime. The softening of the quantum fluctuation spectrum is 
evaluated analytically by the theory of the functional determinants and computed at low T up 
to the crossover. In particular, the negative eigenvalue, causing an imaginary contribution to the 
partition function, is studied in detail by solving the Lame equation which governs the fluctuation 
spectrum. For a heavvy particle mass, the decay rate shows a remarkable temperature dependence 
mainly ascribable to a low lying soft mode and, approaching the crossover, it increases by a factor 



> 

(N 

oo 

f^ . five over the predictions of the zero temperature theory. Just beyond the peak value, the classical 

^^ \ Arrhenius behavior takes over. A similar trend is found studying the quartic metastable potential 



but the lifetime of the latter is longer by a factor ten than in a cubic potential with same parameters. 
Some formal analogies with noise-induced transitions in classically activated metastable systems 
are discussed. 

PACS numbers: 03.65.Sq, 03.75.Lm, 05.30.-d, 31.15.xk 



1. Introduction 

Occuring in a broad class of natural phenomena met ast ability has been intensively stud- 
ied in different research fields ranging from condensed matter to chemical physics, from 
particle and nuclear physics to cosmology. While the mathematics of the metastable phase 
in statistical physics was elaborated by Langer ij] forty years ago, the thorough extension of 
the formalism to the quantum field theory came in the seventies of the last century 2|, y, |4|. 
As the semiclassical approximation is known to be appropriate 5|| to deal with quantum 
tunneling, the path integral method in the Euclidean version [g, |7| has emerged as the fun- 



damental tool for quantitative theories of the decay rate F of a metastable state. The latter 
is given in semiclassical theory by 

r= Aexp{-B/K)[l + 0{K)] (1) 

where A and B depend on the specific shape of the potential. In the path integral 
approach, B is the Euclidean action calculated at the classical path around which the action 
varies the least while A is obtained by summing over the quantum fluctuations paths. In 
quadratic approximation the quantum fluctuations are decoupled and their contribution 
can be evaluated by the theory of the functional determinants which implements the path 
integral method [81]. 

Eq. ([1]) is formally identical to the tunneling energy in the bistable 0^ potential 



which admits an instanton 



id, 



111 ] as the classical path interpolating between the two 
vacua. The fundamental difference lies however in the fact that the instanton minimizes 
the Euclidean action whereas the classical path in the metastable potential, the bounce, 
is a saddle point for the action. Hence the bounce time derivative, which is eigenfunction 
of the quantum fluctuations Sturm-Liouville operator with zero eigenvalue (due to time 
Tanslational invariance), is not the ground state. There must be a negative eigenvalue 



12l ] in the fluctuation spectrum causing an imaginary contribution to the total partition 
function. The decay rate of the false vacuum arises precisely from this negative eigenvalue 
which, moreover, is unique as the bounce path has one node along the time axis. 

While the general formalism of metastability is settled at least in one dimension and 

nn 

in the semiclassical method |13l . I14I |. the evaluations of T for specific systems are generally 



based on the Langer-Coleman approach [15| which assumes a system with infinite size L and 
therefore, strictly speaking, holds in the zero temperature limit of the quantum field theory 
once L is taken equal to Hc/KbT. The same formal approach can be used in classical systems 
to calculate the noise-induced transition rate between locally stable states, for example the 
domain reversal rate in micromagnets [iQ] and the instabilities in metallic nanowires [iTJ. 
In such systems however the character of the transition between activation regimes crucially 
depends on the finite spatial extent [18| , thus quantitative estimates of the Kramer escape 
rate require an extension of the Langer-Coleman method to the more complex finite size 
case. In quantum tunneling systems, the counterpart of the latter is represented by the 
finite T quantum regime below the crossover temperature T* at which thermal activated 
processes set in. Below T*, the decay rate stems from purely quantum fluctuations effects 
and the F dependence on T is governed by the softening of the fluctuation spectrum. To 
determine this phenomenon with accuracy one has to face the non trivial task of computing 



fluctuation functional determinants at low but finite T [l9l] which also includes a knowledge 
of the negative eigenvalue in the spectrum. This paper focuses on this issues developing the 
finite temperature theory of metastability below T* for the case of a model cubic potential. 
In particular, the softening of the fluctuations is analysed quantitatively by means of the 
Jacobi elliptic functions formalism 



20l | which has the advantage to monitor the evolution 
of the system at any T, free from approximations regarding the form of the bounce and its 
fluctuation spectrum. The explicit equations which determine such softening are provided. 



Activation effects 



21 



22| are not discussed in this work which assumes a closed system. 



Metastability in the presence of dissipation has been studied in Ref.[23|. 

As in the space-time Euclidean path integral the time is an inverse temperature, the 
zero T theory assumes that the bounce classical motion lasts an infinite time whereas the 
finite T theory realistically admits that the excursion time for the bounce may be finite. 
This consideration permits to derive the family of energy dependent paths which solve the 
classical equation of motion as shown in Section 2 for the cubic potential. In Section 3, 
using the theory of the functional determinants, I determine analytically and compute the 
quantum fluctuations contribution to the path integral. The solution of the Lame equation 
which governs the fluctuation spectrum is presented together with the evaluation of the 
temperature effects on the lowest lying eigenmodes and their eigenvalues. Section 4 contains 
the results for the decay rate of the metastable cubic potential and a comparison with the 



case of the inetastable quartic potential. Some final remarks are made in Section 5. 

2. Cubic Potential Model 

Consider a particle of mass M in the one dimensional cubic potential: 

V{x) = x^ - -x^ (2) 

plotted in Fig. [TJ^a) for several energies huj. Say a the position of the top of barrier 
whose height is V{a) = 70^/6 with 7 = Mut^/a. I take throughout the paper, a = lA and 
M = lO^rrie with rrie being the electron mass. At a; = the particle is in a local minimum, 
see Fig. [I](a), from which it cannot escape classically and, in the real time formalism, only 
Xci = solves the classical equation of motion. Such local minimum is however metastable 
as quantum fluctuations allow the particle to explore the abyss at x > 3a/2. Let's see how. 

A. Classical Path 

A non trivial classical solution can be found in the Euclidean space. In fact, performing a 
Wick rotation from the real to the imaginary time, t -^ —iT, the equation of motion reads: 



Mx^iir) = V'ix,i{r)) (3) 

where V means derivative with respect to Xd- The particle path x(r) has been split in 
the sum of a classical and a quantum component, x{t) = Xci{t) + ?7(t). The Wick rotation 
is equivalent to turn the potential upside down with respect to the real time as shown in 
Fig. [11(b) for hLj = lOmeV. In the reversed potential the classical motion can take place. 

Integrating Eq. ([3]), one gets: 



^xliT)-Vixdir))= E (4) 

with the constant E representing the classical energy. Defining: 




FIG. 1: (Color online) (a) Cubic potential (in meV) for several oscillator energies, (b) Cubic 
potential for hw = WmeV in the imaginary time representation. The intersections with the 
constant energy E (in meV) define the turning points for the classical motion. 



Xcl[T) = 



_ ^Xciir) 



AE 



27V{a) 



(5) 



Eq. (jl]) is easily integrated to yield: 



To 



±- 



Xci{t) 



dx 



(6) 



where tq is the center of motion between the turning points at which the path velocity 
vanishes. Thus the physical picture set by the boundary conditions is that, see Fig. [11(b), of 
a particle starting from X2 at the time r = —L/2, reaching Xi at r = Tq and returning to 
the initial position at r = L/2. Then, Eq. ([6]) has a time reversal invariant solution whose 
period L is finite and dependent on E. Phase transitions phenomena in spatially extended 
systems are analysed by formally similar models with L being the size of the system 18l. |24|. 



Looking at Fig. [Hb), one captures qualitatively that the amplitude xi — X2 attains the 
largest value for the E = motion while X2 and x^ coincide, x^ becomes negative for E < 
motions. 

The turning points Xi, X2, x^ are given by the zeros of the equation — x^ + X^ + '^ = (x = 
2x/(3a)) which admits three real solutions for k G [—4/27,0], that is for E G [— V^(a),0]. 
After some algebra I find: 



1 2 
Xi= 2 + 2 cos(^) 

X2= i + ^cos(^-27r/3) 

X3 = 3 + 3 cos(^ - 47r/3) 

?9 = - arccos(— - + l) (7) 



At the bounds of the energy range, Eq. (JTj) yields: 



E = ^ Xi = 1; X2 = X3 = 

E = -V(a) ^ xi = X2 = 2/3; xs = -1/3 



(8) 



Thus, at the sphaleron energy Egp^ = \E\ = V{a) |25|, the amphtude of the finite time 



solution has to shrink into a point. Let's find the general solution of Eq. ([6]) by pinning the 



center of motion at Xd{To) = Xi and using the result 



3: 



XI 



dx 2F(A,p) 



:iM Vixi-X){x-X2){x-X3) ^X\ - X3 

X . I ^X\-Xci{j) 
A = arcsm | \ / 

Xi -X2 



P=J^^^ (9) 

V xi - X3 

where F{X,p) is the elliptic integral of the first kind with amplitude A and modulus p. 
Then, through Eqs. ([S]), (jH]), ^, I derive the bounce solution of the finite time theory: 



Xdir) = Y[xicn'^i^,P) + X2sn'^itu,p)] 

.UJ 



(10) 

sn{w,p) and cn{zu,p) are the sine- and cosine- amplitudes respectively [20|. The modulus 
p incorporates the classical mechanics of the problem through the second of Eq. ([5]) and 
Eq. (jTj). At E =0, p = 1, the bounce of the infinite time theory is recovered: 



1.5- 



0.5 



L = 3.37 T =3.44 -E/V(a)=10" 
L = 0.84 T*=13.8 -E/V(a) = 0.1 
L = 0.69 T* = 1 6.8 -E / V(a) = 0.5 



co=10 




1 

X e [ - L / 2, L / 2 ] 



FIG. 2: (Color online) Shape of the bounce solution for three classical energies E. The oscillation 
period L/h is in units meV~^ and T* in K. 



Xci{r) = —cn^{w,l) = ysec/i^(-(r-ro)) 



(11) 



At the sphaleron energy, p = 0, the bounce solution is in fact a point-like object set 
at the bottom of the valley in the reversed potential: Xdir) = a. Thus, Eq. (TTOj) defines 
the transition state which is a saddle for the action below the sphaleron. Representing the 
Jacobi elliptic functions in Eq. (TTUI) through trigonometric series |26| suitable to computation, 
I plot the classical path in Fig. [2] which makes evident how the bounce amplitude contracts 
by increasing the energy over potential height ratio (in absolute value). While the bounce 
can be physically interpreted as a pair of domain walls, it is clear that the extent of the 
r-region over which the domain walls separate the turning points continuously shrinks by 
approaching the crossover. The center of motion tq has been set equal to zero with no loss 
of generality being the system invariant for time translation. 

As the bounce is a combination of squared Jacobi elliptic functions, its period is 2K{p) 
with K(p) = F{tc/2,p) being the complete elliptic integral of the first kind [20(]. Hence, 
from Eq. (ITO|) . I get: 



./j^^^Ys'^L = K{p) (12) 

which estabhshes the relation between the oscillation period and the classical energy 
embedded in the turning points. As in the Matsubara formalism, one can map the imaginary 
time onto the temperature axis, L = H/IKbT*), where T* is the temperature at which the 
particle makes the journey to and from the edge of the abyss for a given E. Then, only 
periodic bounces whose period equals the inverse temperature determine the decay rate and 
the finite time theory can be viewed as a finite T* theory. From Eq. flT2l) I get: 



KsT* = ^^/^ll^ (13) 

The L and T* values given in Fig. [2] are computed on the base of Eqs. (fT2| (fT3|) respec- 
tively. 

At the sphaleron, Eq. flT5]) yields 

huj , , 

KbT: = — 14 

that marks the transition temperature between quantum and activated regimes. Such 
value represents the upper bound for the validity of the model and precisely sets the Goldan- 
skii criterion 27(1 for the cubic potential. Taking hw = lOmeV, the following calculations 



are carried out in the low temperature range up to T* = 18.4697^. 

B. Classical Action 

The classical action A[xc/] for the bounce in the finite temperature theory can be computed 
either in terms of the path velocity Xci{t)'. 



/TO 
dT[x,i{T)f ~ E ■ L{E) 
-L/2 

Q 

Xci{t) = — J^ ■ sn{w,p)cn(zj,p)dn('co,p) 

^ = -'^(Xi - X2) VXi - X3 (15) 



{dn{zu,p) is the delta- amplitude [20|]) or, in terms of the potential V{xc 



A[x,i] = ^^^^ r dxV^ + x'-x' - E ■ L{E) 



(16) 



Both ways require computation of L{E) through Eqs. ([5]), ([7]), (fT2|) . 
In the ii^ — > limit, I get the result 



A[xd] QM^uj 



3, ,5 



(17) 



which serves as testbench for the computational method. The inverse dependence on 
the anharmonic force constant 7 reflects the well known fact that metastable systems are 
intrinsecally non perturbative and provides the fundamental motivation for the semiclassical 
treatment. Eq. (TTTj) permits to set the potential parameters such as the condition A[xci] ^ h 
holds and the semiclassical method is thus justified. As M and a have been taken constant, 
-^yf^ (X CO in the E ^ limit. The bounce velocity and the classical action are displayed 
in Fig. [3] and Fig. H] respectively. Note that A[xci] decreases smoothly with increasing tem- 
perature and absolute value of energy (inset (a) in Fig. jH) suggesting that the transition to 
the activated regime above T* is of second order [l2|, |28| . Although I'm dealing here with a 
closed system, it is worth emphasizing that the inclusion of a dissipative environment 29 1 
would decrease T* but would not change the character of the crossover. 

In general, both first and second order transitions have been found in tunneling systems 



depending on the shape of the potential [30|, ImI, l32|. The criterion used to establish the 



order of the transitions in periodic tunneling systems has been formulated by Chudnovsky 
331] on the base of the behavior of L{E): a monotonic dependence of the period below the 
sphaleron points to a smooth crossover whereas a nonmonotonic L{E) indicates a sharp 
transition. The inset (b) in Fig. H] clearly shows that the former case applies to the cubic 
potential in Eq. ([2]). L{E) decreases monotonically and, consistently, the action versus T* is 
convex upwards. At the sphaleron, I find numerically: L{Esph)/h = 0.628meV~^. Note that 
such value corresponds to 27r/hw 3J] and this is not by accident. In fact at the sphaleron 
the bounce is a point, that is a static solution of Eq. (jl]) but, near the sphaleron, the periodic 
path is the sum of the sphaleron and an oscillation with negative eigenvalue £_! whose period 
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FIG. 3: (Color online) Bounce velocity versus imaginary time for the same parameters as in Fig. [2l 
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FIG. 4: (Color online) Classical action (in units h) versus temperature. Insets: (a) Classical Action 
versus Energy over Barrier potential height ratio on a linear scale; (b) Oscillation period (over h) 
versus Energy over Barrier potential height ratio on a log scale. 



tends to L{Esph) = '^'^/ \/\e-i\ [35]. Then one infers that, for \E\ — > Egph, £-i — > —uj'^. I'll 
come back to this point towards the end of the next Section. 



3. Quantum Fluctuations 



Let's expand the action around the classical path up to second order in the fluctuations 



10 



A[x]^ A[x,i]+5A[r,] + ^5^A[r,] 



MM = / t/r( Mxe/(r)77(r) + V\xa{r))r,{r] 

5'A[ri] ^ 2AM 

^/M= / dT{-f^\r) + -V"{x^,{T))rf{r] 



(18) 



Being the classical solution an extremum for the action, the first functional derivative 
vanishes and, after partial integration in the first term of 5A{ri), one gets the Euler-Lagrange 
Eq. ©. 

Further, differentiating Eq. ([3]) with respect to r, one observes that Xci{t) solves the 
homogeneous equation associated to the second order Schrodinger-like differential operator 
in Af[ri\ 

Or]n{r) = EnVnir) 

6 ^ -dl + V"{x,i{t))/M 

CO 



n= —1 



where the £„ are the fluctuations eigenvalues while the coefficients <j„ of the series expan- 
sion in ortonormal components r7„(r) define the measure of the fluctuation paths integration. 
Then, as a general consequence of the r-translational invariance, Xci{t) is proportional to 
the eigenmode 77o(t) with Eq = 0. The latter cannot be the ground state as, see Fig. [3], Xci{t) 
has one node hence, the quantum fluctuation ground state has negative eigenvalue. From 
Eq. (fT5l) . note that for any two points zui, ZU2 such that W2 = 'cui±2K{p), XcK^a) = idi'^i)- 
The important consequence is that the fluctuation eigenmodes obey periodic boundary con- 
ditions (PBC). 

After these observations we are ready to calculate the space-time Euclidean path integral 
between the positions Xi and Xf connected in the time L. In the semiclassical model, it is 
given by 



11 



< Xf\xi >L= exp 



/OO n Q 
1 J —c 



A[xci] 



h 



Drjexp 



h 



d<in 



'2'Kh/M 



(20) 



K depends only on the functional integral measure. Xi and Xf coincide for the periodic 
bounce, thus Eq. (l20l) represents the single bounce contribution Zi to the total partition 
function Zt- The latter also contains the effects of all multiple (non interacting) excursions 
to and from the abyss which is equivalent to sum [7| over an infinite number of single bounce 
contributions like Zi. Moreover, also the static solution of Eq. ([3]), Xc/ = contributes with 
the harmonic partition function Zh which can be easily determined using the measure in 
Eq. (!20|) . Summing up, Zt is given by 



Zt = Zhexp{Zi/Zh) 
Zh= ^Det[h]\'^/^ 



(21) 



Det[h] being the harmonic fluctuation determinant. Being the decay rate F proportional 
to the imaginary exponential argument through the Feynman-Kac formula [7|], it follows 
that there is no need to determine N as it cancels out in the ratio Zi/Z^. 

Let's proceed to evaluate the fluctuation term in Zi which, formally, is worked out by 
carrying out Gaussian path integrals yielding: 



Drjexp 
Det[d] = 



Aflv] 



h 



n 



^■Det 



O 



-1/2 



(22) 



Eq. (!22|) is however divergent due to the Goldstone mode which reflects the fact that 
To, the center of the bounce, can be located arbitrarily inside L. However such mode can 
be treated separately [281]: the divergent integral over the coordinate d<,Q associated to the 
zero mode in the measure Dr] is transformed into a (Itq integral. Accordingly the eigenvalue 
£q = is extracted from Det[0] and its contribution to Eq. (l22l) is replaced as follows 



12 




(23) 



To be rigorous, this replacement holds in the approximation of quadratic fluctuations pi 
while higher order terms may be significant around the crossover. It is also worth noticing 
that Eq. (12^ is often encountered in the form (eg)"^''^ -^ y 2th -^- However, this is correct 
only in the zero T limit where A[xc/] equals MN~'^ (see Eq. (TT5l) ) while, approaching T*, 
the difference between the two objects gets large. Exactly at T*, the norm of the bounce 
velocity vanishes while the action is finite. This fact is crucial in establishing the behavior 
of the decay rate at the crossover as shown below. 

Now I tackle the problem of the evaluation of the regularized determinant Det^[0] defined 
hy Det[d] = eo-Det^[d]. 

A. Fluctuation Determinant 

The calculation of Det^[0] requires knowledge only of the classical paths. This is a 
fundamental feature of the theory of functional determinants of second order differential 



operators first developed for Dirichlet boundary conditions ^ 



and then extended to general 



operators and boundary conditions in several ways 36|, ISTJ]. As shown above the path 
velocity obeys PBC for any two points tz7i,ti72 separated by the period 2K{p). The latter 
corresponds to the oscillation period L along the r-axis. It can be easily checked that also 
the path acceleration fulfills the PBC. Then, the regularized determinant is given by 



/o(tJ7i)l4/(/o,/i) 

where /o, /i are two independent solutions of the homogeneous equation: 0?7„(r) = 0. 
W{fo,fi) is their Wronskian and < /o|/o > is the squared norm, /q is obviously Xci while 
/i can be taken as 

/.^^;,-/ (25) 
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Using Eq. (fTOl) . derivatives and properties of the elliptic functions, I derive for /i: 



/i 



3a 

y 



■en (ct7,pj + — — sn [w,p) 



dq 



dq 



dxi 
dq 

dx2 
dq 

dFi\,p) 
dq 



-2(xi - X2)sn{zu,p)cn{w,p)dn{ui,p) ■ 



l-p2 



dFi\,p) 
dq 



2(l-p2+p4)3/2 

-1 
2(l-p2+p4)3/2 



1 
2^ 



-E'(A, p) — p F{X, p) sin A cos A 



pz 



V 1 — p2 sin^ A 



-2 1 2 

p = I — p 



(26) 



where E{X,p) is the elliptic integral of the second kind. 

The Wronskian, being constant along r, can be calculated in any convenient point. Let's 
take To as /o(to) = 0. Then: 



_ 9 
~ 8 



-/o(T-o)/l(ro) 



TO 



-a'^uj'^iXi -X2)(xi -Xs) 



dxi 
dq 



(27) 



and < /o|/o >= N'^ given in Eq. (|23l). 

Thus, all ingredients are known to calculate Det^[0]: evaluating Eq. fl26l) on the r— 
range boundaries, Eq. fl2^ transforms into 



Det^[d] 



^a/Xi - XsP 



f^2 



E(7r/2,p)-p2ir(p) 



pz 



< /ol/o > 
W^(/o,/i) 



(28) 



which can be computed using Eqs. (125]) . (127|) . It is however known in the theory of 
functional determinants [8| that only ratios of determinants are meaningful in value and sign, 
such ratios arising naturally in the path integral method as it has been pointed out above. 
In fact, Det^[0] would diverge in the T* ^ limit due to the fact that the determinant 
is the product over an infinite number of eigenvalues with magnitude greater than one. 



14 



Consistently with Eq. (12T]) . Det^[0] has to be normahzed over Det[h] {h = —d'^ + uo'^) 
which, in the case of PBC, is: Det[h] = — 4sinh^(ci;L/2). The normahzation cancels the 
exponential divergence and makes the ratio finite. 
Then, observing that for E ^ ( T* — > 0): 



WifoJ- 



TO o 

6 2 



< /ol/o >^ -a ^ 
5 

K{p) -^ ln(4/v/l-p2) (29) 

from Eq. (!28!l . I finally get the finite ratio 

^^ ^ -^ (30) 

Det[h] 60l^' ^ ' 

The dimensionality [uj~'^] correctly accounts for the fact that one eigenvalue has been 
extracted from Det^[0]. The computation of Eq. fl25]) becomes rather time consuming in 
the T* ^ limit where elliptic integrals need to be evaluated with high accuracy in order 
to get the correct exponential divergence. 

The plot of Det^[0] versus temperature is shown in Fig. [5] up to T* = 18.339/C < T* 
while the inset displays the computed inverse ratio normalized over — GOcj^: the T* ^ limit 
is in fact an excellent estimate up to T* ~ lOK whereas a strong deviation is found at larger 
T* up to ~ T*. Exactly at T*, Det^[0] vanishes thus causing the divergence of the inverse 
ratio. This divergence may look surprising since the zero mode had been extracted from 
the determinant: as it will be unravelled in the next Subsection, there must be a quantum 
fiuctuation mode which softens by increasing T* and ultimately vanishes at the sphaleron. 

1 want to emphasize the key role played by the bounce velocity norm in getting Det^[0] ~ 
at T* ~ T*: close to the sphaleron, < /o|/o > oc p^ and W{fo, /i) oc p^ thus the squared 
norm drops to zero faster than any other term in Eq. ( l28l) . This feature is made evident 
in Fig. E] where the plots of Det^[0]/< /o|/o > and the squared norm are given separately 
against temperature up to T*. 

Eventually, note that the determinants ratio is negative and this sign has physical meaning 
as it is due precisely to the negative ground state of the fiuctuation spectrum. As Zi/Zh 

15 




FIG. 5: (Color online) Regularized fluctuation determinant versus temperature up to T* = 
18.339i^. The inset shows the ratio between harmonic oscillator and regularized fluctuation deter- 
minants. Such ratio is normalized over the characteristic value —QQuj"^ found in the zero temperature 
theory. 





FIG. 6: (Color online) (a) Regularized fluctuation determinant over norm of the zero eigenmode; 
(b) norm of the zero eigenmode. Both curves are plotted versus temperature including the crossover 
point. 

contributes to the partition function by the square root of the fluctuation determinants 
(inverse) ratio it follows that such contribution is purely imaginary. To weigh the effect of 
the negative eigenvalue e^i, its analytical expression has to be determined at finite T* . This 
is done in the next Subsection. 

B. Lame Equation 

Let's take Eq. flT9|) with the second derivative of the potential given by 
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V"{xa{r))= Mu^(i-'^x,i{t)) 



(31) 



Then, using Eq. (TTOl) and working out the algebra, I get the stabihty equation which 
governs the fluctuation spectrum 



(P 



\Vn{T) = [l2p^sn^{zu,p) +An]r]n{r) 
4(1 -3xi) 



dzu' 

A 

Xl -X3 

12 = 1(1 + 1) 



4e. 



^^(Xi - Xa) 



(32) 



This is the Lame equation in the Jacobian form for the case / = 3 [20]. For a given I 
and p, Eq. fl32l) yields periodic solutions (which can be expanded in infinite series) for an 
infinite sequence of characteristic An values. The continuum of the fluctuation spectrum 
stems from this sequence. However, being / positive and integer, the first 2/ + 1 solutions 
of Eq. (1321) are not infinite series but polynomials in the Jacobi elliptic functions with real 
period 2K{p) or AK{p). Being the period of the potential |38|, 2K{p) plays the role of a 
lattice constant. 

Then, Eq. (1321) admits 7 polynomial solutions with eigenvalues An-, n G [— /, /], from which 
the corresponding £„ are derived [39|: 



6^= UJ «! + 02 



e_3 = UJ «! + ^2 



£2=1^ «! + ^2 



5(1 + p^) + 2^Ap^-7p^ + A 
5(l+/)-2v/4p4-7p2 + 4 ^ 

5 + 2p2 + 2^p'^-p^ + A 



£^2= ^ «! + "2 



5 + 2p2 - 2^p^-p'^ + A 



ei= UJ I ai + a2 

,2 



£-1 = UJ ( ai + 0:2 

£0=0 

«! = 1 - 3xi 

Xl ~X3 



2 + 5/ + 2^/Ap^-p'^ + l 
2 + 5/ - 2v/4p4_p2^i \ 



02 



4 



(33) 
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FIG. 7: (Color online) (a) Absolute value of regularized fluctuation determinant over ground state 
eigenvalue ratio (units meV~'^) up to T* = 18.339i^; (b) Ground state fluctuation eigenvalue 
(units meV'^) up to T*. 

The physically relevant feature is that the Sn depend, through p, on the system size L (or 
on T*). However not all the e„ are good fluctuation eigenvalues. In fact the first four values 
have to be discarded as their eigenfunctions do not fulfill the PBC required for the fluctuation 
components: rjnij^i) = Tjnij^i =F 2-ft'(p)). Thus, three good eigenvalues in polynomial form 
are found: i) Eq = is the zero mode eigenvalue which is consistently found also through the 
stability equation, ii) ei lies in the continuum and vanishes at p = 0: this is the soft mode 
responsible for the sharp drop observed in Det^[0] as the sphaleron is approached. At T*, ei 
and £0 merge consistently with the double degeneracy of the corresponding eigenmodes above 
the crossover. Hi) e-i is the negative eigenvalue responsible for the decay of the metastable 
state. In Fig. [71 the ratio Det^[0]/e-i and S-i are plotted separately to emphasize the 
role of the latter. The dependence on T* is obtained as usual by means of Eq. flT51) . While 
the overall shape of Det^[0] is not essentially modified by extracting S-i (compare with 
Fig. [5]), the interesting feature is represented by the softening of e_i with respect to the 
value e_i = — 5u;^/4 found at T* = 0. The substantial reduction (in absolute value) starts 
up at T* ~ lOi^, that is in the same range at which the classical properties deviate from 
the predictions of the zero temperature theory. Finally, at the sphaleron, from Eq. ( l33l) I 
get E-i = —Lo^ consistently with the observation made at the end of the previous Section. 
This completes the analysis of the quantum fluctuation spectrum. 
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4. Decay Rate 

The decay rate in Eq. ([1]) for the finite temperature theory can be now computed as the 



coefficients A and B, identified with \/\Det[h]/Det[0]\/L and A[xci] respectively, have been 
analysed in detail. The general expression for r(T*) is: 



r(r*) = hJ^N'\ 



Detlh] 

exp 



27rh \DetRR[d]\e 



A[x,i] 

h 



(34) 



with Det^^[0] = Det^[0]/e^i. Essentially, the two bound state eigenvalues ^o aiid £_! 
have been extracted from Det[0] and treated separately. 

It is now clear that, in order to calculate the decay rate, there would be no need to 
compute N^^ due to the zero mode J40| as it cancels out with the same term embedded in 
Det^[0] and discussed above. Thus, it is the piece of determinant in Fig. [6](a) which enters 
the computation of F. 

Eq. (!34l) is plotted in Fig. [8] against temperature up to T* = 18.469-/^ together with 
the constant decay rate of the zero temperature theory. There is a substantial increase in F 
above T* ~ lOi^ where the combined effects of quantum ffuctuations and classical action 
softening become evident. Approaching the crossover (at Tp = 18.138i^), F reaches a peak 
value which is larger by a factor 5 with respect to the value of the zero T theory although 
the linewidth remains much smaller than the oscillator energy, T{Tp)/hw ~ 10^^. Certainly, 
such ratio strongly depends on the particle mass through the exponential Gamow factor as 
74[xcz] oc M. Beyond Tp and before T*, the quantum decay rate smoothly merges with the 
classical Arrhenius factor as A[xd]/h — > V{a)/KBT*. In principle this should occur exactly 
at T* once the nonuniform bounce leaves room to the constant solution Xci{t) = a. But 
close to T* the decay rate drops to zero as F(T*) oc (T* — T*)^/"^ signalling that the quantum 



tunneling effects is ending. The power law dependence is driven by \j N~'^ / Det^[0] oc p. 
For this reason the Arrhenius behavior in the decay rate, which is a continuous function of 
T*, takes over slightly before T* and precisely at T\ = 18.25K corresponding to the symbol 
in Fig. El 

The inset shows the plot of the quantum statistical decay rate F^ [l3| which I have 
computed, on the base of the WKB turning point formula 12], by summing over the particle 
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energies in the cubic potential case. In the quantum regime, the calculation converges by 
taking about 1000 energy values (below the sphaleron V{a)) in the Boltzmann average. 
In the activated regime (with energies above the sphaleron), Fs continuously grows with 
exponential behavior while, the matching between the quantum and activated r^ plots 
occurs at around T*, with states of energy E ~ V{a) yielding a significant contribution to 
Fe in the classical regime. 

The deviation with respect to the T = result found in the quantum decay rate F(T*) 
(and culminating in the peak value) is interesting also in view of a comparison with classi- 
cally activated systems in which spatio-temporal noise induces transitions between locally 
stable states of a nonlinear potential. The changes in radius and the stability conditions 



of metastable metallic nanowires are an example of current interest 17|, |4l|, |42|. In these 
systems, the finite size shapes the activation behavior and a power-law divergence in the 
escape rate with critical exponent 1/2 is predicted once a critical lenghtscale is approached. 
In the decay rate of Eq. ( l34|) the small parameter is h which, unlike the noise in classical 
systems, cannot be varied as a function of T (or L). Accordingly the quantum tunneling 
decay rate is always small. Similar conclusions regarding the transition can be drawn by 
carrying out the study of F for the metastable quartic potential 

V{x) = -^x' - -x' (35) 

according to the same semiclassical method described so far. The parameters M, u and a 
have been taken as in the cubic case while 6, in units eVA~^, is given hj 6 = Muj'^/a'^. There 
are in fact some fundamental differences which may be summarized as follows: 1) for the 
quartic potential, the bounce solution is linear in the Jacobi elliptic dn-amplitude and the 
crossover temperature T* is larger by a factor ^/2 with respect to the value given in Eq. (HM . 
2) The classical action is slightly larger than that computed in Fig. H] and, in the E -^ 
limit, v4[a;c/] -^ 4M^co'^/(3M). 3) The harmonic oscillator over regularized fluctuations 
determinants ratio is smaller and, in the T* — > limit: Det[h]/Det^'[0] -^ —VluJ^. 4) The 
Lame equation, which has to be solved for the case / = 2, has the negative eigenvalue scaling 
from E-i = — Scj^ at T* = to E-i = —liJ^ at T* = T*. All these features work together 
to decrease the decay rate of the metastable quartic potential with respect to the cubic 
one. The quantitative results are reported on in Fig. [9] by plotting F for the two potentials. 
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FIG. 8: (Color online) Decay rate (in meV) for a particle in a metastable cubic potential versus 
temperature up to T*. The result of the standard instantonic theory is shown for comparison. 
The symbol marks the point in which the quantum decay rate merges with the classical Arrhenius 
factor. The quantum statistical decay rate Fs is plotted in the inset. 

It appears that, at a given T* within the quantum regime, the life of a particle set in a 
metastable quartic potential is longer (by an order of magnitude) than in a cubic potential 
with the same structural parameters. Close to T*, the decay rate of the quartic potential 
shows the same behavior found in the cubic case as a consequence of the softening of the 
polynomial fluctuation eigenmode with lowest lying positive eigenvalue. 

5. Conclusion 

The semiclassical path integral method has been applied to compute the decay rate of a 
metastable cubic potential at flnite temperature. I have considered a non dissipative system 
where the inverse temperature, introduced as the imaginary time of the Euclidean path 
integral, represents a measure of the system size. Solving the Euler-Lagrange equation for 
the cubic potential I have determined the classical background, the bounce, in terms of Jacobi 
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FIG. 9: (Color online) Decay rates of metastable cubic and quartic potentials versus temperature 
up to T*. 

elliptic functions whose energy dependent modulus keeps track of the classical mechanics. 
The periodicity property of the elliptic solution sets the fundamental relation between bounce 
oscillation period and classical energy. As the former varies monotonically with the latter, 
the transition from the quantum to the activated regime is in fact a temperature driven 
smooth crossover occuring at T* oc uj. Then the quantum regime may cover a broader 
temperature range in systems with characteristic frequencies larger than the value taken 
here. In the latter cases however the quadratic approximation for the quantum fluctuations 
may become questionable. 

In general, the character of this transition in tunneling systems may depend on the form 
of the potential. For the metastable potential the classical background represented by the 
bounce has zero topological charge, fulfills periodic boundary conditions and the amplitude 
of its motion progressively shrinks by increasing T*. At T*, the bounce degenerates into 
a point, the constant solution which determines the classical Arrhenius factor. Then the 
nonuniform bounce samples a T* dependent portion of the (reversed) potential valley and 
continuosly adapts its shape to join the turning points. That's why the crossover is smooth. 

As in the semiclassical method classical and quantum degrees of freedom are intertwinned 
the quantum fluctuation spectrum can be evaluated, by the theory of the functional deter- 
minants, in terms of the classical bounce. More precisely, in terms of its time and modulus 
derivatives. I have obtained a compact expression for the fluctuations contribution to the 
path integral and computed the renormalization of the fluctuation spectrum due to temper- 
ature effects inside the quantum regime. The overall softening is remarkable: it shows up 
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at about T* /2 whereas at lower temperatures the predictions of the zero T (infinite time) 
theory are accurate. The polynomial eigenmodes of the spectrum have been analytically 
determined by solving the stability equation associated to the Sturm-Liouville differential 
operator, the Lame equation. The focus on the ground state is justified by the fact that its 
eigenvalue is negative hence, it is responsible for the decay rate F of the metastable potential. 
The negative eigenvalue softens considerably by approaching the sphaleron where, in turn, 
it determines the oscillation period of the classical bounce. Also the positive polynomial 
eigenvalue plays a key role as it even vanishes at the sphaleron and ultimately induces the 
rapid increase in the inverse determinants ratio. The explicitly derived equations govern- 
ing the temperature dependence of the fiuctuations may be adapted to specific situations 
by tuning the potential parameters. For the specific choice of parameters here presented, 
M = W^rrie and hw = lOmeV, the computed F starts to increase over the zero T decay 
rate at about T*/2 and, close to T*, the peak value is larger by a factor 5, a sizeable effect. 
At the crossover the quantum decay process ends and the Arrhenius classical behavior sets 
in. Finally, a comparison between cubic and quartic metastable potentials having the same 
structural parameters shows that the latter is less metastable in the sense that the lifetime 
is found to be larger by a factor ten at any temperature in the quantum regime. 
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